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Abstract 



We describe a simple and predictive scenario for baryogenesis in theories 
with large extra dimensions which resembles Affleck-Dine baryogenesis. The 
Affleck-Dine field is a complex scalar field carrying a U(l) x charge which is 
dynamically broken after the end of inflation. This generates an excess of 
X over X-, which then decays into Standard Model fermions to produce an 
excess of baryons over anti-baryons. Our model is very constrained because 
the Affleck-Dine field has to be sufficiently flat during inflation. It is also a 
source for density fluctuations which can be tested in the coming satellite and 
balloon experiments. 



I. INTRODUCTION 



Recently it has been proposed that large extra spatial dimensions can explain the ap- 
parent weakness of the electroweak scale as compared to gravity in 3 + 1 dimensions. In 
such a scenario four dimensional world is assumed to be a flat hypersurface, called a brane, 
which is embedded in a higher dimensional space-time, which is known as the bulk. The 
hierarchy problem is then resolved by assuming that TeV scale can be the fundamental scale 
in higher dimensions [0-0]. This however requires the size of the extra dimensions to be 
large. Such a large volume can substantiate the hierarchy in energy scales. The volume 
suppression Va, the effective four dimensional Planck mass M p , and, the fundamental scale 
in 4 + d dimensions M* are all related to each other by a simple mathematical relation 

Ml = Ml +d V d . (1) 
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This automatically sets the present common size of all the extra dimensions at bo. For two 
extra dimensions, and, M* = 1 TeV, the required size is of order 0.2 mm right on the current 
experimental limit for the search of deviations in Newton's gravity 01]. Recent experimental 
bounds suggest M* to be much larger. Naturally, such model has an important impact 
on collider experiments ||, and on cosmology. In this paper we address one of the most 
important issues in particle cosmology, the origin of baryon asymmetry in theories with large 
extra dimensions. 

The generation of baryon asymmetry requires three well-known conditions; C and CP 
violation, B or L violation, and out of equilibrium decay ||. It is quite probable that the 
early Universe had strong departure from thermal equilibrium due to large expansion rate in 
the early Universe. However, achieving out of equilibrium condition becomes more difficult 
when the inflationary scale and the final reheat temperature of the Universe is lowered down 
to the electroweak scale and below. As we shall see, this is a major obstacle for realizing 
baryogenesis in the context of large extra dimensions. Another major problem is a proton 
decay. A low fundamental scale induces fast proton decay via dimension 6 baryon number 
violating operator in the Standard Model (SM). With a low fundamental scale the usual 
coupling suppression is not sufficient. 

The cosmological setup for large extra dimensions is quite different from the conventional 
one. Firstly, if the electroweak scale is the fundamental scale in higher dimensions then there 
can be no massive fields beyond the electroweak scale in four dimensions. Secondly, the size 
of the extra dimensions can be quite large as compared to the electroweak scale, which 
implies a new degree of freedom with a small mass scale related to the size of the extra 
dimension. This field is usually known as the radion field. Its mass can be as small as 
0(eV) for two large extra dimensions. The stabilization of such large extra dimensions is a 
dynamical issue because they grow from their natural scale of compactification ~ (TeV) -1 
to the millimeter size in order to solve the hierarchy problem. In fact, the stabilization 
must take place at the very initial phase of the Universe via some trapping mechanism as 
discussed in Ref. 0. 

Another challenge is how to realize inflation in these models such that one naturally 
obtains the correct amplitude for the density perturbations. There have been many proposals 
||, but, the most appealing one is invoking a SM singlet scalar living in the bulk ||, and we 
shall point out why this is the only mechanism which works well. There is also the question 
of the presence of Kaluza Klein (KK) states of the graviton and any other fields residing 
in the bulk. At high temperatures it is quite possible to excite these KK states. Above 
a certain temperature known as the normalcy temperature, the Universe could be filled by 
the KK modes. The Big Bang Nucleosynthesis (BBN) constrains the normalcy temperature 
to be above ~ 1 MeV. In order to be safe from other cosmological bounds the final reheat 
temperature which is the largest temperature one can envisage during the radiation era, 
should be smaller than the normalcy temperature, which is constrained by cosmological 
considerations to be as small as 100 MeV 10 . 



There are many ways of generating baryon asymmetry in the Universe, one of which is the 



simplest and predictive scheme known as the Affleck-Dine mechanism [11], first discussed in 



the context of supersymmetry. A scalar condensate which carries non-zero baryonic, or/and 
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leptonic charge survives during inflation and decays into SM fermions to provide a net baryon 
asymmetry. The adaptation and elaboration of this particular mechanism in theories with 
large extra dimensions is the main goal of this paper. It has already been shown that the 
Affleck-Dine mechanism is the only solution for providing an adequate baryon asymmetry 
in this context with a very low reheat temperature, see Ref. . 

We begin our discussion by introducing some of the salient cosmological features of 
the large extra dimensions. In section II, we briefly discuss the consequences of having a 
temperature larger than the normalcy temperature. This is necessary in order to judge the 
merit of other baryogenesis scenarios such as leptogenesis and electroweak baryogenesis. In 
section III we discuss inflation and its couplings to the SM fields. We shall also discuss 
very briefly the radion stabilization because of its importance. In section IV we consider 
various other proposals of baryogenesis and argue that within the context of large extra 
dimensions, where the inflaton is a higher dimensional field, they essentially fail to provide 
the observed baryon asymmetry. In section V we introduce our Affleck-Dine baryogenesis 
model which we embed within an inflationary scenario. In section VI we describe various 
facets of our model and argue why it appears to be the sole mechanism which can generate 
the correct baryon asymmetry. We shall also highlight various ways of constraining the 
model parameters. Note even though in most of the cases we specialize for the critical two 
extra dimensions, our applications and conclusions are valid for any number of extra spatial 
dimensions. Finally, in section VII we summarize our main results. 



II. THERMAL HISTORY OF THE UNIVERSE 

In order to appreciate thermal history of the Universe we need to estimate the tem- 
perature below which the Universe could be safely regarded as a radiation dominated one, 
which we call here the normalcy temperature, T c . In this section we shall review some of 
the already known results and provide some new insight on how the cosmological evolution 
changes if the temperature of the Universe exceeds T c . It is possible to excite the KK states 
from the plasma with a temperature T. The production of these states from relativistic 
particles depends on the cross section which is Planck mass suppressed. It is noticeable that 
the cross section of processes such as 7 + 7 — > G is of order ct 7+7 ^g ~ (TR) d /M^ [j], where 
R is the effective size of the extra dimensions. Once gravitons are produced their evolution 
can be traced by the Boltzmann equation 

— + 3Hn G = (a v)^ ~ ^ - — , (2) 

where tq is the decay life time of a massive KK mode, given by 

M? 
m 



^3- ( 3 ) 



If we naively assume that right after the end of inflation the first term is dominating the 
right-hand side of Eq. then by taking a(t)T(t) = constant, where a is the scale factor 
of the Universe, we can in fact simplify the above equation. While doing so we may also 
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neglect the evolution of the individual mode and shall concentrate upon all possible KK 
states excited up to a given temperature. We obtain 

d(n G /n^) {a v)n 7 

dT ~ HT ' U 

Once the KK states are excited they are no longer in thermal equilibrium. However, indi- 
vidual KK modes with a mass ~ T remain present in the Universe until they decay. We 
can integrate Eq. (|4]) while assuming that we are in a standard cosmological era such that 
H 2 oc p/M p , to get 

n G (T) n 7 (T r )(q v) 
n 7 H(T r ) ' { ' 

The temperature T r designates the largest temperature during radiation era, known as the 
reheat temperature of the Universe. We also take v = 1, henceforth. Now substituting the 
cross section and assuming that the relativistic particles dominate the Universe; n 7 ~ T r 3 , 
we evaluate the right-hand side of Eq. (j5|). The ratio thus obtained can not exceed more 
than one at any later times in order to maintain the successes of nucleosynthesis era and so 
we obtain a simple bound on T r , which is given by [l[ 

(M d+2 \ l,l+d 

Note, for the preferred values M* ~ 30 TeV for d — 2, we obtain T c < 100 MeV. This 
is an extremely strong constraint on the thermal history of the Universe. It reiterates 
very strongly that the radiation dominated Universe simply can not prevail beyond this 
temperature. Therefore, any physical phenomena such as first order phase transition, out- 
of-equilibrium decay of heavy particles, if at all taking place beyond T c , need to be revised. 
Now, it is pertinent to ask how the evolution of the Universe changes if we assume that it 
is possible to exceed the normalcy temperature, but let us first note that as we take larger 
number of extra dimensions, d — ► oo, we apparently increase T c up to the fundamental scale 
T c - M* 0. 

A priori it can not be ruled out that the temperature of the Universe should not exceed T c . 
The reheat temperature is determined by the inflaton coupling to the matter fields. The only 
constraint upon T r is that it must be larger than few MeV. Therefore, there is no guarantee 
that the Universe can not have a radiation domination with an instantaneous temperature 
more than the normalcy temperature. In such a case, there is a plenty of KK states which 
can be excited easily beyond T c , which actually leads to enhancing the number of relativistic 
degrees of freedom. The number of relativistic degrees of freedom increases as g(T) ~ (RT) d , 
which also determines the number of degrees of freedom determining entropy. At this point 
one may wonder how the KK modes, which actually have Planck mass suppressed couplings, 
can be brought into thermal equilibrium. Indeed, they cannot if n G <C n^. However, we 
are in a opposite limit when the KK gravitons have started dominating the number density 
of relativistic decay products of inflaton. Once the KK modes are produced they go out 
of equilibrium due to the fact that the self interaction among these gravitons is extremely 
weak. However, at the time they are being produced they introduce an extra entropy to the 
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already existing plasma. The entropy injection becomes important once we are above T c 
(the KK modes are excited with masses up to T > T c ). The distribution function for each 
and every KK mode has a completely different profile, for the KK mode with mass m ~ T c 
can be understood to be close to the relativistic particles, but the same can not be true for 
the heavier KK modes. As a result the effective entropy stored within the KK modes can 
be written as 

s = ^ ~ (RT) d T 3 , (7) 

neglecting factors 0(1). The term in the bracket corresponds to the relativistic degrees of 
freedom which is roughly the number of degrees of KK modes. We should mention here 
that while deriving Eq. (^), we have naively assumed that the final distribution function for 
the KK modes are peaked around the final temperature and all the KK modes below that 
are produced abundantly with an uniform distribution which mimics that of a relativistic 
species. 

Now, following the fact that entropy conservation gives 

sa 3 = T d+3 a 3 = constant , (8) 

we obtain a simple relationship between the expansion of the Universe and the rate of change 
of the temperature 

H = -* + **, (9) 
3 T K ' 

which actually determines the cooling rate of the Universe provided the KK modes are fairly 
stable. This then leads to an approximate Hubble parameter 

T (d+A)/2 r>d/2 

H(T) ps — . (10) 

However, the KK modes do decay and some of them would perhaps decay much before 
nucleosynthesis, but in our case this depends on the maximum temperature one could reach 
beyond the normalcy temperature. Note, in a standard cosmology with a radiation dom- 
inated Universe the temperature scales like inverse of the scale factor. In such a case any 
particle which has a Planck mass suppressed couplings with other fields can decay before 
nucleosynthesis provided its mass is beyond 10 TeV. Otherwise, lighter particles decay much 
after Nucleosynthesis and their number density is well constrained from the diffusion of 



gamma rays in a micro wave background radiation ||13|| . However, if we had a KK domi- 
nated Universe, then there might be a temperature-scale factor relationship which would be 
governed by Eq. (|]). This would also affect the life time of the massive KK modes. The 
simplest estimation for the masses of the KK modes which would decay before nucleosynthe- 
sis could be obtained by demanding that a KK mode with a mass just above the normalcy 
temperature should decay before nucleosynthesis. This determines the temperature of the 
Universe at the time of the decay. For instance, for d = 2 extra dimensions we obtain 

Tdecay OC R' 1 , (11) 
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where R is the size of the extra dimensions. In order to be consistent, the above temperature 
ought to be more than ~ 0(MeV). This restricts the size of the extra dimensions to be much 
smaller than ~ 10~ 10 mm. This translates into a lower bound for the fundamental scale which 
is now increased to M* > 10 8 GeV. This result is inconsistent with our basic assumptions 
of having a low quantum gravity scale. This certainly may rule out the possibility of having 
large extra dimensions for d = 2. The situation improves a little if there are more than 
two extra dimensions. For instance, if we take d = 6 the above analysis gives a bound for 
the size R < 2 x 10~ 14 mm, and, for the fundamental scale M* ^2x 10 5 GeV. Since, the 
normalcy temperature also increases, which can be ~ 0(1) TeV, baryogenesis may not be 
so troublesome provided we also take into account at least six extra dimensions. However, 
a simple hope like this seems to be a mirage once we realize that inflaton must reside in 
the bulk whose couplings automatically determine a reheat temperature below electro weak 
scale. This conclusion is actually quite robust and regardless the number of extra dimensions. 

Our result asserts that we must reheat our Universe below T c given by Eq. in order 
not to excite the KK modes with an over abundance. We have noticed that normalcy 
temperature is the largest temperature of the Universe below which one can safely regard 
the content of the Universe as a radiation dominated one. 



III. INFLATION, REHEATING AND THERMALIZATION 

A. Inflation and density perturbations 

In order to provide a relatively small reheat temperature one automatically requires very 
small couplings of the inflaton to the matter fields. As we already know that the dynamics 
of the inflaton plays a crucial role in achieving out-of-equilibrium condition for baryogenesis. 
Keeping this in mind we briefly review inflationary dynamics which shall also act as a preview 
for the AD baryogenesis. 

As discussed above, the largest temperature in a radiation dominated Universe must 
be smaller than 100 MeV for two large extra dimensions. One way to achieve this is to 
assume that the inflaton is living in the higher dimensional bulk, and upon compactification 
it naturally admits Planck suppressed couplings to the SM fields. A simple model of inflation 
can be constructed using only coupled scalar fields in 4+d dimensions || (for other attempts, 
see Ref. ||). The potential can be written down as 

V(N, 0) = \*Mt (nZ - ^A 2 ) 2 + ^ + , (12) 

where <f) is the inflaton field, and, N is the subsidiary field which is responsible for the 
phase transition. The coupling constants are g and A, in general they are different, and Nq 
determines the vacuum expectation value. Note that the higher dimensional field has a mass 
of dimension 1 + d/2, which leads to non-renormalizable interaction terms. However, the 
suppression is given by the fundamental scale instead of the four dimensional Planck mass. 
Upon dimensional reduction the effective four dimensional fields, <p, N are related to their 
higher dimensional relatives by a simple scaling 
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N 



'V A N. 



(13) 



From the point of view of four dimensions the extra dimensions are assumed to be compact- 
ified on a d dimensional Ricci flat manifold with a radii b(t), which have a minimum at bo . 
The higher dimensional metric then reads 



d s 2 = dt 2 - a 2 (t)dx 2 - b 2 (t)df 



(14) 



where x denotes three spatial dimensions, and y collectively denote the extra dimensions. 
The scale factor of a four dimensional space-time is denoted by a(t). After dimensional 
reduction the effective four dimensional action reads 



S 



J d A X\f^g 



Ml 



1 



-R + -d.ad^a -U(a, , ,,,, 
16vr 2 2 



0, 



+ -d.Nd^N 
2 ^ 



-exp(-da/a )V((f>,N)) 
The potential V(<f>, N) can be derived from Eqs. (|T2D and (|T3|). 



<P 2 N 2 + X -m\<? 



(15) 



(16) 



The radion field a(t) can be written in terms of the radii of the extra dimensions 



a(t) = cr In 



m 

b 



00 



d(d + 2)M 2 



16n 



1/2 



(17) 



From the above equation, it is evident that o"o is proportional to the four dimensional Planck 
mass. For b(t) ~ (TeV)^ 1 , and, bo ~ 1mm, the modulus of the radion field takes a very 
large initial value. The radion field has a potential which at the minimum is given by 
U(a) ~ m 2 ,a 2 , where m T ~ !CT 2 eV. In this paper we shall assume that such a potential is 



essential to stabilize the large extra dimensions (for discussion, see Ref. [ft!]) 



There are mainly two phases of subsequent inflation in this scenario. We notice that 
initially N is stuck in the false vacuum determined by 



V 



A 2 iV 4 + -m 2 ^ 2 



(18) 



X 2 M 2 M 2 if N 



M m , the contribution is almost constant if the 



which renders V(4>, N) ,~ 
vacuum dominates over the second term in Eq. (|T8|), which is necessary for the generation 
of density perturbations. The first phase of inflation occurs when <fr is still rolling down the 
potential and has not undergone a second order phase transition. During the first phase of 
inflation the exponential term due to the radion field present in front of the scalar potential, 
see Eq. fll5]) is responsible for driving a power law inflation. For details we refer 0. The 
exponential potential inflates not only the brane but also the bulk, which expands from 
(TeV)" 1 to mm size, which is equivalent to 35 e-foldings of inflation if there are only two 
extra dimensions. Towards the end of this phase when the radion field a — ► ctq, the radion 



gets a running mass m 2 s 



0(1)H 2 + m 2 . 



The Hubble parameter H 2 3> m 2 , 



therefore 
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the Hubble correction dominates the radion mass. Once, the radion gains this mass it 
simply rolls down to its minimum as dictated by U(a). This ensures that not only the 
extra dimensions have already been grown to the adequate size 60 ~mm, but also have been 
stabilized. This takes place while the inflaton is still rolling down the potential and N is 
still locked in the false vacuum. From the last term in Eq. ( pT5| ) one finds that the exponential 
exp(cr/(Xo) becomes of 0(1) when a is trapped in its minimum, and the effective potential 
is solely given by Eq. (TTS]). Therefore the dynamics of a is completely frozen except for 
quantum fluctuations. 

In a completely model independent way it can be easily argued that at least 43 e-foldings 
of inflation is necessary in order to produce structure formation . The amount of required 
inflation is less than the usual 60 e-folding because of a smaller inflationary scale and reheat 
temperature. The scale of inflation during this second phase is simply given by 

A.Yj 




(19) 

If N ~ M*, we automatically get the Hubble expansion at the end of inflation: H « M*. 
The adiabatic density perturbations can be generated while the inflaton is rolling down 
the potential until it reaches a critical value where a second order phase transition takes 
place and the effective mass square for N field becomes negative. The critical point can be 



found from the potential Eq. (|i8|) to be 

. _ A (N M p 



r g V M* 

Note, if A ~ g, and, Nq ~ M*, we automatically get C ~ M p . 

The two slow-roll conditions for inflation can be laid down very easily 



(20) 



2 



M p 




(21) 



16tt \ V((f) } N) 

M= 87 V&N) = 8^livfJ {\N ) <<X ' (22) 

where prime denotes derivative with respect to 0. The spectrum of the density perturbations 
is given by |L5| 

32V^JV) ± (23) 
H 75 Ml e 43 ' { ' 

where 643 is defined roughly 43 e-foldings before the end of second phase of inflation. From 
data <5h — 1-91 x 10" 5 The interesting cosmological scale leaves the horizon when 

043 = c e 43 '' 7 ' ~ 0c, provided we assume \r)\ 1/43. By proper substitution we get a simple 
relationship 



s 



In order to illustrate, let us set Nq = M* ~ 100 TeV and ~ 10 GeV; we then obtain 

5h~(A 2 (7)x10- 5 , (25) 



which provides the right COBE (and Boomerang) amplitude with A ~ g ~ Oil). Let us 
stress that all the above conclusions are independent of the number of extra dimensions. 

Let us mention here that we actually require Nq < M* for two reasons. First, N provides 
a mass squared contribution to the Higgs potential, which has to be less than the electro- 
weak scale in order to provide the right magnitude for the Higgs mass. The second obvious 
point is that effective mass for and N must be greater than the Hubble rate H after the 
end of inflation. This is necessary to terminate inflation just after the phase transition. 
However, for most of our calculation Nq ~ M* remains a very good approximation in order 
to show the merit of the baryogenesis model. 

The spectral index n = l + 2r] — 6e is presently constrained by observations \n — 1| < 0.13. 
If we naively assume Nq = M#, and, A = g, we find 

AM* ~ 1.27 x 10 15 M GeV . (26) 

This suggests that \r)\ has to be extremely small. Given that, e is also very small, see Eq. (pT|). 
we conclude that our model predicts a perfect scale invariant density perturbations. If the 
coupling is of order one and M* ~ 100 TeV, the slow roll parameter 

\v\ <10~ 10 . (27) 

If we assume that the inflaton sector is solely responsible for the adiabatic density pertur- 
bations, the small slow roll parameter constrains the amplitude of any other scalar fields we 
intend to introduce in our setup. This we shall elaborate when we discuss AD baryogenesis. 



B. Post inflationary dynamics of (j) and N 

Let us first discuss the classical dynamics of the fields <p and N after the phase transition. 
During this era both the fields oscillate with an initial amplitude determined by 

The frequencies are determined by the effective mass scales of the fields at the global mini- 
mum, given by 

m = V2gN , m N = V2XN . (29) 

We notice that the phase transition must terminate inflation immediately, otherwise there 
could be another bout of inflation which might give rise a particular signature in the spectrum 
of the microwave background radiation which we will discuss in a separate publication. 
A slow transition might provide density perturbations of order one and might produce 
primordial black hole formation. In this paper we do not take into account of this phase. 
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Therefore, in order to ensure ~ > H 0} the constraint on N < yo/87rM*, for A ~ g. 
However, in order to illustrate our model for baryogenesis the essential physics remains for 
N ^M*. 

When the frequencies of the oscillations are different the fields loose their coherence and 
the motion becomes chaotic, see Ref. |T7|-[nj. If the coupling strengths g and A are equal 
we can easily obtain the field trajectories around the bottom of the potential, which follows 



a straight line in a phase space of — iV [|18ipJ| . There exists a particular solution of the 



classical equations of motion given by [0,0 

N(t) = V2(<j )c -<j ) (t)) . (30) 

It can be argued that in a static Universe it is possible to obtain a classical solution for either 
or iV near the bottom of the potential. Initial motion is an-harmonic, but the oscillations 
near N Q can be expressed in a simple form: 



N(t)^V2(Mz)N 



MJ-"l 1+1 3 COS{m * t] 
In an expanding Universe the amplitudes of the oscillations decay and 



— — cos(m0 t) 



(31) 



(32) 



where the amplitude of the oscillations decreases as $(£) oc 1/t. The precise form of $(£) 
depends on the ratio H/rh^. If the ratio is large, the decay of the amplitude is felt in a 
couple of oscillations, otherwise, it may take many oscillations before the expansion leads 
to decaying in amplitude. In our case we can make a rough estimation. We notice that 
~ mjv ~ H at the time of phase transition. Therefore, the fields begin to roll down 
towards their respective minimum of the potential from their initial amplitude which is 
~ M p . The amplitude of the fields decreases very quickly until it reaches a point when 
fh^ ~ ffiN > H . This happens because the Hubble parameter H is also decreasing as the 
fields roll down as ~ 1/t, where t is the physical time starting from the end of inflation. In 
order to see this, we need to ensure 

Hit) _ 1 



which happens when 



< 1 , (33) 



1 2vriV osc . , 

— = — (34) 

where N osc is an approximate number of oscillations. For ~ H , we obtain N osc > l/2ir. 
Therefore, within one oscillation the Hubble parameter decreases quite rapidly. This makes 
sure that the fields having masses fh$ ~ fh^ ~ Hq can oscillate about their minimum with 
a common frequency for many oscillations before they completely decay. We may estimate 
the number of oscillations to be iV osc ~ fh^r = (M?/M%) ~ 10 28 . Hence during the initial 
stages the dominant effect of expansion will be to render the total oscillations more harmonic 



around the minimum. Therefore, our assumption of the field evolution given by Eq. (32) 
holds well. However, in reality we need the opposite limit on masses, which would only 
slow down the decay of the amplitude of the oscillations. With this brief discussion on the 
dynamics of the fields we move on to discussing thermalization. 
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C. Reheating and thermalization 



It is believed that the total energy density of the inflaton is transferred into radiation. 
The minimal requirement is to have a thermal bath with a temperature more than 0(1) MeV 
in order to preserve the successes of BBN. Recall that now we can not reheat the Universe 
above normalcy temperature T c . The final reheat temperature depends on the decay rates 
of the oscillating field || 



T r ~o.i v /(r ( ^ + r w )M p . (35) 

Therefore, we need extremely weak, non-renormalizable coupling of and iV to other SM 
fields. In our model this is natural because the bulk fields have Planck suppressed interac- 
tions with matter fields stuck on the brane. However, this fact also causes some problems, 
e.g. why the zero mode inflaton is reheating the brane and why not the bulk? The infla- 
ton could decay into some other lighter degrees of freedom into the bulk. This point has 
been already addressed in Ref. || for the case of gravitons, and here we recapitulate the 
arguments. 

In our previous section we found that both the fields get an effective mass term ~ 
rh N oc v^-XiVo, if we assume A = g. Setting N ~ C(TeV), then both and N are 
kinematically allowed to decay into Higgs field, h, with a mass (9(100GeV). This is different 
from the usual Kaluza-Klein theories where the production of matter through inflaton decay 
occurs everywhere. The decay rate is estimated as follows f| 

f 2 M* 



32vrM2m . 



N 



where, / is the coupling constant that we take of order one. If rh^ ~ fn^ ~ 0.1M*, then the 
reheat temperature is T r < 100 MeV, which is more than 0(1) MeV and below the normalcy 
temperature, even for the most constrained case of d = 2 extra dimensions. However, note 
that such a low reheat temperatures is a generic prediction if the inflaton field is living in 
the bulk. 

The next point is concerning the production of KK graviton via the inflaton decay. The 
KK gravitons can be directly produced from the KK modes of and N, and possibly the 
KK modes of other scalar fields present in the bulk via n — > <piG n _i interactions, where 
n, I are the KK numbers, and G is the KK graviton. On the other hand n modes can be 
produced via collision processes, such as 00 — > n 0_ n , and similar reaction for N. The rate 
of exciting the KK modes and their subsequent decay rate to graviton can be estimated as 
1 

F m n mj . . 

where rr? n = + n 2 /R 2 is the excited KK mode. Note that the iViV or 00 scattering rates 
for producing their KK counterparts are smaller than the direct decay of 0, iV to the brane 
fields. Therefore, the zero modes of and iV still prefers the Higgs as a final decay product. 
However, late during the thermalization era the scattering phenomena may give rise to the 
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production of KK modes of inflaton. In the above equation we have only estimated a single 
process for graviton production. There are plenty of other accessible modes in the final 
channel. This enhances the decay rate 

IV, total = IV-xfcGn-1 ~ ,n L 2 • ( 38 ) 
I * 

The excited KK inflaton or KK partners of other scalar fields present in the spectrum are 
extremely short lived. The heavier KK mode decays into the lighter KK mode plus gravitons, 
and eventually all the KK modes of the inflaton decay into the zero mode. In case where 
extra bulk fields are present the reheating of the bulk can be naturally avoided, if either 
those modes are as heavy as the inflaton, or, the effective inflaton couplings to the bulk 
fields are smaller than the inflaton-brane interactions. 

So far we have tacitly assumed that reheating is almost instantaneous. This might not be 
the case, especially when the fields are oscillating. During this period, the equation of state of 
the Universe in most of the cases is given by that of a matter dominated era. However, if 0, N 
are decaying very slowly, then they might also decay into lighter degrees of freedom, such as 
relativistic species directly. One might expect that this could be the most preferred channel. 
However, this is not the case, because again the oscillating fields have non-renormalizable 
couplings to these lighter fields and the decay rate follows r^,jv-^ 77 ~ M^/M^. Notice, its 
resemblance with that of Eq. (|36|), if fh^ ~ M*. This is an important lesson and all it tells 
us that 0, N decaying into Higgses and into lighter degrees of freedom is equally preferred. 
However, if there were some radiation which could thermalize, then in such a case, the 
Universe could in principle follow an equation of state which would be determined mainly 
by a mixture of relativistic species, inflaton, and non-relativistic Higgses. This would only 
affect the thermal history of the Universe, and not the dynamics, because inflaton energy 
density is dominating over all. It has been shown in the simplest situation where inflaton and 
radiation components are allowed, that the instantaneous temperature of a plasma might 
exceed the reheat temperature of the Universe. In order this to happen one must also satisfy 



H >T (f, + T n fll3fl . We notice that this condition may be satisfied in our case at the very 
beginning of the oscillations, but in spite of large initial amplitudes for 0, N, the oscillations 
are damped quickly. 

The temperature of the plasma may reach its maximum when a/a ~ 1.48, where a 
denotes the scale factor of the Universe and the subscript denotes the era when inflation 
comes to an end. The maximum temperature is given by |20 



^ a77 (^J W) • m 

where g* denotes the effective number of relativistic degrees of freedom. For the purpose 
of illustration if we fix Hq ~ M* and T r ~ 100 MeV we found T max ~ 10 5 GeV. This 
temperature is much higher than the actual reheat temperature, and it seems that this is 
a generic prediction of inflationary scenarios in extra dimensions. However, this warrants a 
preferred production of relativistic species from the decay of inflaton. This is unfortunately 
certainly not the case with the present situation. In order to proceed with our present 
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discussion, we note that after reaching the maximum temperature, the temperature of the 
intermediate plasma decreases as 

T~1.3^^^y /4 T max a- 3 / 8 , (40) 

The important thing to notice here is that the thermal history of the Universe is again 
different: the temperature does not drop like T oc cT 1 so that the entropy of the plasma 
s oc a 15 / 8 . This will eventually dilute any KK graviton being produced during this era, and, 
as long as ng/n 7 < 1, we would not expect any further alteration of the results mentioned 
in Eqs. fl3|) and (0). 



IV. BARYOGENESIS 

The constraints on inflationary parameters which we have been discussing so far must 
be bear in mind when discussing baryogenesis. We certainly need a concrete mechanism 
from particle physics in order to address baryogenesis. The problem here is that we should 
introduce some baryon number violating processes. However, in theories where the funda- 
mental scale is low naturally introduces dangerous higher order operators. For instance, 
dimension 6 operators such as QQQL/A 2 , where Q's and L correspond quark and lepton 
SM doublets, which can mediate proton decay, unless A > 10 15 GeV. This operator violates 
baryon and lepton number conservation by AB = AL = 1. There is also a possibility of 
having right handed singlets uuude/A 2 , for which A > 10 12 GeV. There are other processes 
which violate baryon number, such as neutron anti-neutron oscillations with a dimension 7 
operator uddudd / A 3 , which implies A > 10 5 GeV. An alternative is the dimension 9 operator 
QQHQQH/A 5 with A > 10 4 . (The experimental bound on proton life time is r p > 10 33 
years |2l|], and for neutron anti-neutron oscillations T nfi > 1.2 x 10 8 seconds [22. 23] . ) 



We assume that the above mentioned operators can be avoided in some way or other. It is 
very difficult to come up with a model where baryon violating operators are not constrained 
at a TeV scale. Thus, one has to ensure that such operators are not being reintroduced 
by the mechanism of baryogenesis. Especially, in our case we ought to be careful with an 
operator such as xQQQL- In case x develops a vacuum expectation value ~ M*, fast proton 
decay is inevitable. 

Moreover, in the context of extra dimensions, leptogenesis is not a viable mechanism as 
we shall argue now. In leptogenesis a net lepton number is produced in the decay process of a 
heavy fermionic singlet such as a right handed neutrino, which is then processed into baryon 



number by anomalous B + L violating sphaleron interactions |24j] . However, the electro weak 



sphaleron transitions are active only up to 100 GeV |25]. In our case viability of leptogenesis 
has a simple catch. A singlet right handed neutrino can naturally couple to the SM lepton 
doublet and the Higgs field through yLHN . This leads to a potentially large Dirac mass 
term unless the Yukawa coupling y ~ 10~ 12 . Moreover, now the see-saw mechanism fails 
to work, since, the largest Majorana mass we may expect can never be larger than the 
fundamental scale. Therefore, given a very small neutrino mass ~ y 2 (H) 2 /M ~ y 2 ■ 0(1) 
GeV, we still have to fine tune y 2 < 10~ 10 , in order to obtain the right order of magnitude 
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for the neutrino mass. In any case the decay rate of the right handed neutrino is very 
suppressed and it is similar to Eq. (|36|). This means that when the right handed neutrino 
decays into the SM fields, the background temperature is of order of the reheat temperature 
~ (9(1 — 10) MeV. At this temperature the sphaleron rate is exponentially suppressed, which 
is actually a set back for leptogenesis. A reheat temperature of at least 0(1 — 100) GeV 



is required for making this scenario viable p6[ , which could in principle be attained if the 



number of extra spatial dimensions is increased to six, at least from the point of view of the 
normalcy temperature. However, as we observed before, in the class of models we considering 
this is not the case. Indeed, as mentioned above already, unless we increase the fundamental 
scale, the largest reheating temperature we can get from the inflaton decay is just barely 
about 100 MeV, regardless the number of extra dimensions. This makes leptogenesis even 
more difficult. 

A different possibility that sphalerons can reprocess a pre-existing charge asymmetry 
into baryon asymmetry |S7| reflected in an excess of over anti-e^ created during inflaton 



oscillations. This mechanism requires (B + L) violating processes to be out of equilibrium 
before cr comes into chemical equilibrium, such that the created baryon asymmetry could 
be preserved. Again, this has to happen at or above 100 GeV. One could then assume 
that inflaton decays preferably into relativistic species such that the plasma thermalizes to 
a temperature above 100 GeV p8| , |29[| . However, this scenario cannot be implemented in the 



context of large extra dimensions because in this case the oscillating inflaton field injects 
more entropy to the thermal bath as discussed earlier. This leads to an immediate dilution 
of whatever baryon asymmetry has been created prior to the reheating era. The dilution 
factor is given by : 

-i = ( s(T t ) \ = ( g.(T r ) \ ^ \ 3 (Mm 3 ^ 

where s is the entropy and Tew ~ 100 GeV. For a lower reheat temperature such as T r ~ 1 
MeV, the above expression gives 7 _1 ^ 10 25 (assuming T oc a" 3 / 8 and g*(T c ) « g*(T T )). 
Therefore, one concludes that the initial rib/s has to be extremely large > 10 15 , this is 
certainly an extraordinary requirement on any natural model of baryogenesis and practically 
impossible to achieve. Even if we increase reheat temperature T r ~ 100 MeV, we would still 
require n^/s of order one. 

We have learned two important lessons. First, the large entropy production during the 
last stages of reheating can in principle wash away any baryon asymmetry produced before 
electroweak scale. Second, it is extremely unlikely that leptogenesis would work. The only 
simple choice left is to produce baryon asymmetry directly. The sole mechanism which can 
be successful in these circumstances seems to be the Affleck-Dine (AD) baryogenesis [IT 
which we are going to discuss now. 



V. AFFLECK-DINE BARYOGENESIS 



For the details of the AD baryogenesis we refer the readers to Ref. 
relevant questions are following: 



30 . In our case the 
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(1) Can we have a condensate ? 

Unlike in the usual AD mechanism, which is based on Minimal Supersymmetric Standard 
Model, we do not have flat directions automatically in-built in our model, or, protected by 
super symmetry. Even if we invoke some flat directions, it cannot be associated with a 
condensate carrying B or L. The condensate can not be protected alone due to lack of any 
symmetry argument. One can not form a SM condensate in the bulk because bulk gravity is 
color and flavor blind. This suggests that we necessarily have to assume some fundamental 
scalar field. 

(2) What charge should it carry ? 

The AD field, which we denote here by x has to be a gauge singlet carrying some global 
charge under U(l) x . This global charge has to be broken dynamically in order to provide 
a small asymmetry in the current density. This shall be reflected by generating an excess 
of x over X- The charge associated with this AD field must be such that baryon number is 
violated maximally. All that we need is to ensure that the SM quarks maintain the small 
asymmetry between baryons and anti baryons. Notice that x field as suc h does not create 
the baryon asymmetry. The asymmetry is produced due to the difference in number density 
of x over X- This is the most important aspect of our model. This small asymmetry is then 
transferred via the decay of x an d X into the SM quarks and leptons. This constraints the 
decay channel for x an d X which we discuss next. 

(3) What kind of interactions should it have ? 

This is a non-trivial issue because we do not have a condensate made up of SM quarks and 
leptons. We need to assume that x interactions with SM fields conserve U(l) x symmetry. 
Therefore, the quarks and leptons must carry a non zero global x charge. However, we do 
not want the x~X asymmetry to be transferred to non-baryon number violating interactions 
such as interactions involving the Higgses. Therefore, the Higgs field should not carry a 
global x charge, forbidding x decay into Higgses. 

Regarding the decay channels, coupling of x to SM dimension 3 operators such as QQ, 
which has AB = 0, cannot provide baryon asymmetry. Similarly, for higher dimensional 
operators such as Qhq, where Q (q) is the right (left) handed SM quark and h is the Higgs 
doublet. A dimension 5 operator cannot be constructed at a quark level because of the 
color symmetry. The lowest order turns out to be the dimension 6 operator QQQl, for 
which AB = AL ^ 0, which can certainly transmute any asymmetry in x to the quark 
sector. Thus, the global charge carried by chi has to be chosen such that QQQl carries the 
opposite one. Thus, forbidding the presence of this operator alone on the theory and making 
the coupling xQQQl the lowest possible order for x- Also note this operator can mediates 
proton decay too, unless one ensures that x does not develop a vacuum expectation value, 
which is a severe constraint on model building, but not a difficult one to realize. In the same 
spirit one may check those operators which induce n — n oscillations. Again, an effective 
AB = 2 operator UDDUDD, or, such as (QQQh) 2 cannot be induced by, say, integrating 
out x- To avoid the decay of x into Higgs fields we take h to be chargeless under U(l) x . 

(4) Can AD field be either <f>, or, N ? 

Given the constraints, neither N nor the inflaton field can act as an AD field. The 
auxiliary field N develops a vacuum expectation value after the end of inflation, which would 
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immediately induce proton decay. Regarding 0, the bad news is that; since the asymmetry in 
AD field depends on the initial configuration, a large amplitude oscillations in can induce 
an undesirably large — asymmetry. Hence we invoke a separate AD field, a fundamental 
scalar field x w hh a global charge U(l) x . 

We remind the readers that the inflaton energy density p\ should govern the evolution 
of the Universe. Eventually the decay products of the inflaton should be responsible for 
reheating the Universe. The inflaton decays before x decays via baryon violating interaction 
and generates a baryon asymmetry given by 

r^r^p, (42) 
s n x m x pi 

The final entropy released by the inflaton decay is given by s xs p\/T r . The ratio rib/n x 
depends on the total phase accumulated by the AD field during its helical motion in the 
background of an oscillating inflaton field. In our calculation we shall always approximate 
the total phase ~ 0(1). 

In theories with extra dimensions there are two choices for the AD field. It could either 
be a brane field or a bulk field like and N. If we assume that the AD field is a brane field, 
then it cannot have an effective mass higher than the fundamental scale, and the effective 
initial field amplitude \x\ < M*. These constraints suggest that the energy stored in x can at 
most be p x ~ m^M^, where m x < M*. The energy density stored in and N is quite large 
because these fields can have initial amplitudes close to Planck scale so that p\ ss M^M^. 
Thus we would find 

— ~ — ^ » 10- 34 ^ « 10- 10 . (43) 
s M p M p M* K J 

for T r ~ (9(1 — 10) MeV, we conclude that x cannot be a brane field. 

Let us therefore promote AD field to the bulk whence the energy density stored in the 
AD field is p x ~ m x M^ since now x can have large initial amplitude |x| > M*. This implies 
that the maximum baryon to entropy ratio is given by 

1~ K (m) (m») ~ 10_1 ° (l&sv) ' «-» 

where we assumed T r ~ 10 MeV and ~ 10 TeV for concreteness. Thus we can achieve an 
adequate baryon to entropy ratio from x decaying into SM quarks and leptons. Of course, 
the numeric value of the ratio depends on the initial amplitude for x an d its mass. Note, if 
we set m x ~ M*, we obtain the initial amplitude x ~ M GUT ~ 10 16 GeV. 

Now we can also estimate the decay rate for x field and x fields. The interaction is given 

by 

xQQQL (45) 



M*M P 



We sum over all possible channels; such as various color and family index combinations 
which can be of the order of thousand. On the other hand we strictly assume that inflaton 
is decaying into Higgses. Hence we estimate the decay rate as 
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r ^(?) 2 (s) ,r — (46) 

In order that the decay products of X thermalizes before Nucleosynthesis we have to assume 
m x pa M*. This does not leave much freedom for the masses which makes the model more 
predictive but also demands some level of fine tuning if AD baryogenesis is to be successful, 
as we shall discuss below. 



VI. THE MODEL 

Let us now describe our model. We assume that the inflaton sector is responsible for 
breaking U(l) x charge dynamically. The AD potential can be written as 

V/ AD (0, N, Xl , X2 ) = k\ 2 N\ X l + xl) + ^ (^) 2 (xl + xl? 

+4(^)\n( X 1- X 1), (47) 

where /«i,/«2,«3 are constants, and Xi an d X2 are the real and imaginary components of 
the complex field X . Note, all the terms are Planck mass suppressed, because Eq. (|47|) 
depicts an effective four dimensional potential derived from higher dimensional Lagrangian 
by integrating out the extra spatial dimensions. From Eq. (f4~7|) , it is evident that since 
during inflation the auxiliary field is trapped in the false vacuum N = 0, this renders Xi an d 
X2 massless and the AD potential becomes almost flat. We also notice that the last term 
is also not present during inflation. Therefore, the only contributions comes from the self 
coupling which obviously allows the AD field to evolve. The details of the field evolution 
shall be discussed in the coming subsections. Note that we do not include a term of the 
type (M* / Mp) 2 <fi 2 \x\ 2 as it would ruin both inflation and baryogenesis because of the large- 
amplitude oscillation of X ; this is a fine tuning we have to perform on our model for which 
we have no dynamical explanation. 



A. Constraining the initial amplitude for the AD field 

Dynamical evolution of multi scalar fields induces density perturbations of two kinds; 
adiabatic and isocurvature. Since, in our case xi an d Xi are massless during inflation they 
contribute to the adiabatic fluctuations. In order to obtain a scale invariant spectrum, the 
fluctuations should be mainly generated by and N. However, this imposes an upper bound 
on the allowed amplitude for the AD field at the time when the modes relevant for large 
scale structure formation are leaving the horizon during inflation, which happens in our case 
during the last 43 e-foldings of inflation. 

In the regime when the main inflationary potential is due to the inflaton sector, the 
slow-roll parameter \t]\ can be written as ]31] 



87r V((f), N) 8tt V(4>,N)V'(4>,N) ' 1 ; 
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where prime in Vad denotes derivative with respect to AD field x- According to Eq. 
the second term in the above equation must be smaller than 10~ 10 . This leads to 

| X (0)| 5 < 327rxlO- 10 x-j(^J M p 5 . (49) 

If A, K2 ~ Oil), and iVo ~ M*, we obtain an initial amplitude for |x(0)| < 10~ 2 M P . This 
is a constraint which we should bear in mind while estimating the total baryon asymmetry. 
Note, for a reasonably fast phase transition Nq < M*, the amplitude for x(0) becomes even 
smaller than the above limit. 

Similarly one might expect isocurvature perturbations arising from the AD field. Notice 
that the angular direction of AD field is effectively massless, but gains a dynamical mass 
> H just after inflation. This is unlike the case of a supersymmetric AD baryogenesis where 
the field remains massless even after the end of inflation provided there is no supergravity 
correction to the U(l) violating term. Nevertheless, in our case we would expect isocurvature 
fluctuations during inflation because there are more than one dynamical scalar fields present 
during inflationary phase. It is well known that in the case of isocurvature fluctuations it is 
easier to constrain the spectral index J}2 ] , rather than the amplitude of perturbations. This 



restricts the initial amplitude for the AD field very much in a similar vein as in Eq. (|49|). 
However, since we are assuming that the main contribution to the density perturbations 
is coming from the adiabatic sector, rather than the isocurvature fluctuations, we would 
not expect any further improvement on the limit we have already obtained upon the initial 
amplitude for xi an d Xi from Eq. 



B. Dynamics of AD field 

The dynamics of x is complicated. Even though x has no effective mass during inflation, 
it has an effective quartic self coupling which determines its dynamics (see Eq. (f4~7D ). The 
equations of motion in terms of the component fields Xii Xi are given by 

Xi + 3H X i = -4(^j (xi+xl)Xu (50) 

X2 + SH X 2 = -4 (^j (xi + xl)X2- (51) 

We remind that Xi an d X2 are the components of a a bulk field, therefore, they would have 
naturally taken initial values close to the fundamental scale in higher dimensions but close 
to the Planck scale in four dimensions. From the point of view of four dimensions, the fields 
simply roll down from the Planck scale because the self coupling induces a curvature terms 
for xi an d X2 fields. We note that the suppression in the couplings is very small ~ (M*/M p ) 2 
if K2 ~ 0(1). The fields are very weakly self coupled. However, in the process of rolling 
down the potential their effective running mass becomes of ~ H, given by Eq. (p^9|). When 
this happens their dynamics is effectively frozen at a particular amplitude which observes 
the constraint derived independently in the earlier section, see Eq. ([49|). 
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As we have discussed, we may assume that xi an d X2 collectively follow a similar trajec- 
tory. The dynamics of |x| freezes when 3iif|x| pa — /^{M* / M p ) 2 \x\ 3 ■ This immediately gives 
us a simple solution for the largest values of the fields: 

where N e = HAt is the total number of e-foldings, which could be larger than ~ 130. The 
total number also takes into account of inflation during the radion stabilization, but for 
the structure formation it is the last 43 e-foldings of inflation which matters. For a simple 
estimate, if A ~ k 2 , we obtain |x(0)| < 10 18 GeV, in well accordance with Eq. fli9|) . When 
inflation comes to an end all the fields begin oscillations; for and N, the initial amplitudes 
are given by Eq. fl2~8D, while for xi and X2 the initial amplitudes are determined by Eq. (f|9j). 



The post-inflationary dynamics of Xi an d X2 depends on <ft and N, which are both 
oscillating. As a result xi an d X2 become massive fields (see Eq. ((4"7|)). This leads to the 
equations of motion 

Xi + SHxi = -2/6? 2 N\, - 4 ^ 2 (xl + xl)Xt ~ 24 (^j 2 0iV Xl , (53) 

X, + 3H X2 = -24 (^j 2 N\ 2 - 4 (^j 2 (xl + X\)X2 + 24 (^j 2 0iV X2 . (54) 

There is an effective mass term for x fields which is again field dependent. The last term in 
the above equations comes with an opposite sign and this is responsible for giving rise to a 
small asymmetry in x over X- As a simplest approximation, we neglect the back reaction of 
X fields on the background fields, which is actually true as long as N, <fi are responsible for 
generating the final entropy of the Universe. This assumption simplifies the situation and 
allows us to estimate x asymmetry created by the motion of Xi an d X2 analytically. 

We have earlier pointed out that <fi and iV oscillate several times before they decay at 
the reheating temperature. While the initial oscillations are quite large ~ M p , the energy 
density stored in the oscillations is quite small due to the presence of very small couplings of 
order (M*/M p ) 2 ~ 1CT 28 . Likewise, Xi an d X2 a lso oscillate with a large amplitude and they 
also go through several oscillation periods before they decay with a rate given by Eq. ([16]). 
By substituting Eq. (|2~8|) in Eq. (f47|), we may estimate the effective masses 



m 



v^/tiiVo , m X2 ss v^/tiiVo . (55) 



The coupling constants K\, K2 can be tuned to obtain m Xl ~ m X2 pa in order to match 
the decay rate of T x ~ r^, ~ T^. This finally ensures that the AD field decays along with 
and N fields, and this also prevents the decay of and iV into x fields. If we compare 
Eqs. (P5|), and fl55]), we notice that all the masses are of same order with a small variation 
due to different couplings. If we set K\ > n 2 ~ A ~ g, we can ensure that xi an d X2 
go through more oscillation periods than the background fields cf) and N. Here we tacitly 
assume that the oscillations in x do not have any dynamical back reaction on cf) and N 
oscillations. This is true if the field couplings are weak and indeed, in our case the couplings 
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such as k 2 (M*/M p ) 2 , k|(M*/M p ) 2 are extremely small even if we choose k± ~ K3 ~ Oil). 
Such a weak coupling also ensures that the field trajectories do not give rise to a chaotic 
behavior. In order to avoid any chaotic behavior it is necessary that all the fields oscillate 
with a similar frequency around the bottom of their respective potential. Hence we choose 
Ki = A = g while is constrained by the expressions Eqs. (f^), and (|52"D, and K3 shall be 
constrained by the baryon to entropy ratio. However, in order to have a helical motion for 
Xi and X2, we ought to have Ki > k 3 , and k 2 > k 3 . 

We still have to estimate the classical behavior of Xi an d X2, which do not develop any 
vacuum expectation value. Their classical evolution can be estimated by 

Xi(t) ~ Xi(0)A(t) cos(m Xl t) , X 2(t) ~ X2^)A{t) cos(m X2 t) , (56) 

where the amplitudes Xi(0) an d X2(0) are constrained by Eqs. ( fl9|) and (|52"D. Both the 
amplitudes decrease in time as A(t) oc 1/t. The frequency of the oscillations are determined 
by Eq. (|5^). In what follows we shall always assume that the relative phase between the 
fields Xi an d X2 is of order 0(1). We now have all the tools to address the generation of 
X ~ X charge asymmetry from the classical evolution of all the fields. This shall provide us 
with a final baryon to entropy ratio in this particular model. 



C. Baryon to entropy ratio 

The CP violation in our model is given by the last term in Eq. ([17]). For a charged scalar 
field this is equivalent to C violation. The B violation arise via the decay of Xi because 
the decay products have AB ^ 0, and we have a non-trivial helical oscillations in x which 
accumulates net CP phase which is transformed into asymmetric x- The net x asymmetry; 
n x can be calculated by evaluating the Boltzmann equation: 

n x = \ (x*X - X*X) ■ (57) 
With the help of Eqs. (|53|) and (|54|) , we can rewrite the above expression as 

n x + 3Hn x = 4«jj ( ^ (iV(t)0(t)) Xl (t) X2 (t) . (58) 

The right hand side of the above equation is the source term which generates a net x 
asymmetry through a non-trivial motion of Xi an d X2 fields. We integrate Eq. ( |58| ) from t 
which corresponds to the end of inflation up to a finite time interval. 

n x a 3 = 44 (^j J t y(t')mW(t')Xi(t')x2(t')dt' . (59) 

The upper limit of integration signifies the end of reheating. We assume that the oscillations 
continue until the fields decay completely. Before we perform the integration, we notice that 
the integrand decreases in time. This can be seen as follows ; first of all notice that the 
approximate number of oscillations are quite large before the fields decay ~ (Mp/M*) 2 . This 
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allows us to average and N oscillations; from Eqs. (|30|) and (|3^) we get the time dependence 
{4>N) ~ $ 2 (t) ~ 1/t 2 . Similarly, from Eq. (p6|) , (X1X2) ~ l/^ 2 provides another suppression. 
While taking care of the expansion, where the scale factor behaves like a(t) ~ t 2//3 , the overall 
behavior of the integrand follows ~ 1/t 2 . This suggests that the maximum contribution to 
X asymmetry comes only at the initial times when to ~ 1 /Ho, where Hq is determined by 
Eq. dl9|). The right hand side of the above equation turns out to be 



"x* 27*3 Ho ■ m 

We have assumed that the total CP phase, which is given by two factors: an initial phase 
determined arbitrarily during inflation and the final dynamical phase which is accumulated 
during the oscillations, is of the order ~ 0(1). We have also assumed N ~ M*. 

The final ratio of \ number density produced and the entropy is given by 



s 27A 2 \ M p J 

where we have used the fact that s oc a 3 in our case. This is the final expression for 
X asymmetry produced during the helical oscillations of \ and is to be compared with the 



observed baryon asymmetry, whose range is 4(3) x 10 <t\b = tib/ti- < 7(10) x 10 [13 



The upper bound on the baryon to entropy ratio given by Eq. (|6l"l) depends on the amplitude 
of Xi an d X2, which were frozen during inflation, it also depends on the number of e-foldings, 
which could be ~ 100. In order to evaluate Eq. (|6~l"D , we may take an example: T r ~ 100 MeV, 
and M* ~ 100 TeV. This gives an asymmetry of order ~ 10~ 10 , if we choose the couplings 
of order one. The final asymmetry in x is transferred to the SM quarks via baryon number 
violating interaction mentioned in Eq. P5|). The asymmetry is injected into thermal bath 
along with the inflaton decay products. It is essential that the thermalization takes place 
after AD field has decayed in order not to wash away the total baryon asymmetry. 

Finally, we mention that our model is quite generic. The robustness of the model is that it 
can work for arbitrary number of extra spatial dimensions. The model predicts baryogenesis 
just above nucleosynthesis scale. The model does not rely on first order phase transition of 
electro weak baryogenesis, nor does it depend on sphaleron transition rate. Our approach 
is similar to Ref. |33|], where this model has been embedded in a supersymmetric theory in 



four dimensions. Like in the case of AD baryogenesis ||31|| , our setup can be tested by the 
forthcoming accurate measurements of the spectral index of the microwave background; our 
model predicts a very flat spectrum. 



VII. CONCLUSION 

In summary, we have presented a natural mechanism of baryogenesis in the context of 
low quantum gravity scale with large compact extra dimensions. Our mechanism is generic 
and it is independent of the fundamental scale and the number of compact extra dimensions. 
This mechanism does not rely on any extra assumption other than invoking a fundamental 
scalar field that lives in the 4 + d dimensional space time. The baryogenesis scheme can work 
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at any temperature lower than the electro weak scale. This is because of the presence of non- 
renormalizable couplings which automatically reheats the Universe with a temperature close 
to BBN. This feature is independent of the number of extra dimensions. Our mechanism 
does not depend on the rate of sphaleron transition and also does not rely on leptogenesis. 

The three requirements for baryogenesis is fulfilled in our case as follows. The C and 
CP violation is a dynamical process in our case. The AD field has a charge associated with 
a global U(l) x symmetry. This symmetry is broken once the AD field gains a dynamical 
mass from its coupling to the inflaton sector. This happens once inflation comes to an end 
and when all the fields start oscillating around their respective minima. The dynamics of 
the AD field is solely responsible for CP violation and the initial CP phase is determined 
during inflation which we have assumed to be of order one. During post inflationary era 
the real and the imaginary parts of the AD field has different equations of motion and as 
a result the motion of AD field is helical in the field space. During this helical motion the 
AD field accumulates a net asymmetry in \ over X- This asymmetry is then transferred to 
the SM sector via a baryon violating interaction which does not allow fast proton decay. 

The decay of inflaton into Higgses and their subsequent decay generates a net entropy 
in the Universe. The thermalization of the Universe is quite late and the final reheat tem- 
perature is as low as 100 MeV for a fundamental scale as small as 100 TeV. The AD field 
decays into quarks and leptons directly imparting the total baryon asymmetry it has gen- 
erated during the helical motion, obviously, the net CP phase generated during this helical 
motion is understood to be also of order one. We note that the baryon to entropy ratio is 
determined by the ratio of the reheat temperature and the fundamental scale. The ratio also 
depends on the number of e-foldings. This tells us that the evolution of the massless AD 
field during the inflationary stage also counts in producing a net baryon asymmetry, through 
fixing the initial amplitude of the AD field. The most important aspect of our model is that 
the model parameters are very constrained from the density perturbations during inflation. 
The initial amplitude of the AD field is also constrained in order not to depart from the scale 
invariance of density perturbations measured by the satellite and balloon experiments. The 
constraint which we apply on the amplitude of the AD field comes purely from adiabatic 
perturbations. However, presence of many fields automatically introduces a possibility of 
isocurvature fluctuations which we shall work out in future. In this sense, our model can be 
falsifiable. As a final note we mention that our dynamical mechanism of producing baryon 
asymmetry can be used even in four dimensions and the mechanism is capable of generating 
baryon asymmetry at a scale much lower than the electro weak scale. 
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